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CALCULATIONS  OF  STEADY  AND  UNSTEADY  TRANSONIC  FLOWS  OVER 


* 


AIRFOILS  BY  MONOTONE  SWITCH 

Ang  Haisong  and  Cheng  Zhonglu 

Shengyang  Research  Institute  of  Aerodynamics 

1.  INTRODUCTION 

When  calculating  transonic  unsteady  potential  flows, 
scientists  are  continuously  seeking  improved  methods  for 
solving  non-linear  problems  which  include  shock-wave 
movement,  and  for  increasing  the  efficiency  of  calculation. 
In  1977,  Ballhaus  and  Goorjian  [1]  succeeded  in  solving  the 
time-field  unsteady  transonic  small-disturbance  equations  by 
using  alternative-direction  implicit-difference  method 
(ADI),  and  this  method  was  called  LTRAN2.  Reference  [5]  also 
carried  out  a  study,  which  is  similar  to  LTRAN2,  by 
introducing  a  coordinate  transformation.  Since  then,  several 
authors  [2]  have  confirmed  through  numerical  experiments 
that  the  suitable  frequency  range  of  the  LTRAN2  method  can 
be  expanded  by  using  a  time-correlation  term  to  improve  the 
original  quasi-steady  boundary  conditions.  In 
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1981,  Goorjian  and  Buskirk  [3]  improved  the  Murman-Cole-type 
switch  by  using  the  Engquist-Osher-type  monotone  switch,  and 
proved  that  its  efficiency  is  higher  and  the  stability  is 
better.  They  eliminated  the  possibility  of  existing 
expansive  shock-wave  non-isentropic  solutions  which  were 
permitted  by  the  original  switch. 

In  order  to  increase  the  efficiency  of  calculation  and 
expand  the  applicable  range,  we  have  introduced  an  E-O-type 
ADI-di fference  switch,  the  body-surface  and  tail-vortex- 
surface  boundary  conditions  with  an  additional  time- 
derivative  term,  and  a  pressure-coefficient  formula.  In 
order  to  reduce  the  region  of  calculation  and  avoid  the 
inconvenience  of  coordinate  transformation,  the  unsteady 
non-ref lective  boundary  conditions  [4]  and  steady  asymptotic 
solution  were  used.  A  set  of  programs  for  calculating  steady 
and  unsteady  flows  using  main  equations  having  the  same 
forms  were  established.  The  numerical  experiments  have  shown 
that  this  is  a  cost-effective  method  for  calculating  two- 
dimensional  transonic  unsteady  flows  with  certain  precision 
and  applicable  scope. 
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2.  BASIC  EQUATION  AND  BOUNDARY  CONDITIONS 


The  main  equation. 

The  equation  of  two-dimensional  unsteady  transonic 
smal l-disturbance  flows  (under  the  assumption  of  low 
frequency)  can  be  written  in  the  following  form: 

(-2kMamV.),+  [(  1  -MoL)*.- r  +  1  +(*,),-  0  (  1  ) 


Here  the  coordinates  x,  y,  time  t,  and  the  disturbance 

speed-potential  are  already  non-dimensionalized  by  c,  o_1 

and  cU*  (c  is  the  chord  lengthi  u  is  the  frequency),  and  the 

simplified  frequency  k=oc/U«.  Exponent  m  is  taken  according 

to  the  series  expansion  analysis  of  the  precise  formula  of 

critical  pressure  coefficient  of  isentropic  flows,  that  is: 

2Y  +  1 

the  corrected  value  which  satisfies:  l  >  *  > 

The  boundary  conditions  on  the  body  surface.  (>* /(*,  t  )) 
The  condition  stating  that  the  relative  velocity  of  air 
flow  is  tangent  to  the  instantaneous  body  surface  is: 


-f.+kf, 


For  the  resonant  situation,  the  boundary  condition  at 
the  wing  surface  can  be  written  as: 

*,-h.-at+ a,«ii(mO+«,*(  *  (  3  ) 
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where  h=h(x)  is  the  geometric  shape  of  body  surface ,  «o 
is  the  attack  angle  at  the  average  position,  ctj  is  the 
amplitude  of  unsteady  movement. 

The  pressure  coefficient  of  unsteady  flows  is: 

C,«-  2(T.+*V,)  (4) 

The  boundary  conditions  on  the  surface  of  the  tail- 
vortex  : 

According  to  the  Kutta  condition,  on  the  surface  of 
tail  vortex,  we  have: 

0  (5) 

AC,-r.+*r,-o  (6) 

where  the  circulation  of  the  tail  vortex  |~(x  ,  t )  =A*(x ,  t )  tai  i  • 
Formula  (6)  reflects  the  characteristic  of  the  downward 
propagation  of  the  tail  vortex. 

The  far-field  condition  of  unsteady  flows: 

At  the  far-upstream  and  far-downstream  positions,  we 
use  the  Engquist-Ma jda  non-ref lective  far-field  boundary 
condition 

1 1  —Mai—i  Y  +  1  )#<£*.!*,-  0 
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At  the  upstream-front,  we  take  *=0  at  the  positions 
which  are  far  enough  and  go  to  the  boundary  of  the  upstream- 


front  by  interpolating.  Downstream  ,  under  the 
assumption  that  the  disturbance  pressure  is  zero,  we  have: 


0 


(8) 


The  equation  and  conditions  for  steady  flows: 

The  main  equation  is  the  same  as  equation  (1).  Only  the 
first  term  coefficient  of  the  iterated  parameter  is 
changed.  The  boundary  conditions  at  the  body  surface  and  the 
tail  vortex  surface  are  those  given  above  by  eliminating  the 
derivative  term  of  t.  The  far-field  condition  is  an 
asymptotic  solution  of  the  velocity  potential  of  the 
disturbance . 

)  +  t»n  1  (■j^')]+~«jr  **+0V  J *  <9> 


where  «=  ^  \  •  and  F(*s)  is  the  thickness 

function  of  the  body  surface. 

The  initial  condition: 

For  the  initial  condition  of  unsteady  flows  we  take  the 
result  of  steady  flow  calculation: 


*<*,  >,  0  )-*>«*(*,>) 


(10) 
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3.  DIFFERENCE  SWITCH 


After  the  differencing  of  equation  (1),  the  equation 
can  be  rewritten  as  an  alternate  implicit  form  which  uses 
double-step  calculation  for  each  time  step: 

In  x  direction: 

2hM  ai(  -  DJ  +  (  1 

In  y  direction': 

2  (  1 

Here : 


A/-3./+S./  (12a) 

~f.  _i_  1.Kl/2(*M-*lH)r  (12b) 

■  |  #  1  "  j  »  I 

A./-(/.  s  )C  1 /2  )3‘*  (12c) 

.  1  "  *  ’  1  1  ~  i  •  * 

1/2  Cc'u  ._l_  . +  (c.+c,a  X*i.i-KhX */-*<.,)*')  (  1 2d  ) 

i  1  j  •  1  •  j  •  / 


1  ^a“*  ci““(  T  +  l  )A/oC,  5"  — c,/c,  (  Sound  Speed  of  disturbance)  (12e) 


J 

*  i 

I*. 


1  .»  S)  ( ^-n  &x?) 

i .  _i_  . !  1  "T*  1  _N 

'  "  I  *  1  I  “min  («  .  i  . •  “)  . 

'  ‘  *  *  1  (in  Axf ) 


(  1 2  f  ) 


jS 

¥ 

V 


h.V.i-W.i-Val.x.,)C  1  /  2 


(  1 2g ) 


( 12h ) 


*  I  .  *6 W.J 


( 12i  ) 
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Here,  *ix  and  *ly  are  first  order  backward  difference 
operators.  The  above  E-O-type  correlation  switch  has  the 
following  monotonic  character: 


(*»rv«w  ,1.)>  0 

•  SI 


i 

•S 


therefore  it  has  better  stability  and  convergence. 

The  switch  difference  of  the  non-ref lective  far-field 
boundary  is: 


! 


!_L  , . .  y i  *  '* • 


$ 

1 


■  >’.•  ■  ■  ' 


»  '  .  0  r  V 


By  using  different  "following  methods"  to  obtain  the 
solutions,  the  four-  and  three-diagonal  algebra  equation 
groups  can  be  obtained. 

4.  EXAMPLES  AND  ANALYSIS 

Steady-flow  calculation: 

We  have  calculated  the  steady  flows  for  several  wing 
types  and  several  Ma«.  Fig.  1  shows  the  calculated  results 
and  a  comparison  for  one  of  the  calculations.  It  can  be 
seen  that  the  results  of  this  method  is  closer  to  the 
experimental  data.  The  standard  for  convergence  is:  the 
maximum  difference  of  velocity  potential  10-4-;  the 

iteration  takes  300  ~  520  steps.  The  non-uniform 
network  of  the  calculated  space  is  (44X35)  and  (50X35)  in  a 
physically  rectangular  coordinate  system,  I  ?  c ,  |  x  L.,^5  c  . 


,iW~ 


“C*'V*;'y v?  •'*"*:»?;« 


(Fig.  2  shows  the  situation  at  a  certain  moment,  and  it  is 
similar  for  other  moments.)  Therefore  we  have  found  that 
with  or  without  the  term  of  the  derivative  of  t,  not  only  is 
the  range  of  frequency  correlated,  but  also  the  correctness 
of  the  calculation  will  be  influenced. 


Fig.  2.  The  pressure  coefficient  of  a^NACA64A006-type  wing 
with  1/4  vibrating  side  wings. 

1  -  this  method;  2  -  this  method;  3  -  boundary  conditions 
contain  the  term  of  the  derivative  of  t;  4  -  convergence 
solution  when  the  boundary  condition  does  not  include  the 
term  of  the  derivative  of  t. 
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Fig.  3.  The  unsteady  pressure  coefficient  of  a  NACA64A006- 
type  wing  with  vibrating  side  wings  (the  shock-wave  aoveaent 
is  B-type) . 

1  -  this  aethod;  2  -  reference  [1];  3  -  shock-wave 
position;  4  -  no  shock  wave;  5  -  no  shock  wave;  6  -  side¬ 
wing  position. 


LTRAN2  has  siaulated  three  types  of  the  periodic 
aoveaent  shock  waves  found  by  the  experiaent  of  Tijdeaan 
[6].  But  this  is  only  a  qualitative  siaulation  (see  Table 


1  1 


1). 


Table  1. 


Experimental  and  calculated  conditions  of  oscillat; 
of  wing  flap  with  NACA64AQ0b  airfoil. 


J 

A 

B 

1/ 

c 

•  * 

Mb 

* 

L  -- 

*  1 

* 

Tifkw  «•* 

II 

•  Ul 

•  IT* 

•  1*4 

•  u 

•  :«• 

LTRANz  tt«$ 

•  m 

•  «M 

•  •'« 

•  IU 

I  til 

•  «M 

1  -  type  of  shock  wave  ■oveient;  2  -  the  Tijdeman 
experiment;  3  -  1TRAN2  calculation;  4  -  A-type;  5  -  B-type; 

6  -  C-type. 

In  the  simulation  of  the  three  shock-wave  movement  with 
this  method,  we  have  found  that  when  the  Ma.  and  k 
conditions  are  directly  taken  from  the  Tijdeman  experiment  , 
The  same  curve  pattern  of  shock-wave  movement  resul ted  ( Fig . 
3  gives  the  results  for  B-type),  and  the  results  are 
compared  with  the  experimental  results  (Fig.  4).  Therefore, 
the  simulation  of  this  method  is  consistent  quantitatively 
with  the  experimental  results. 


Fig.  4.  A  companion  of  the  calculated  results  and  the 
expenaentsl  results  of  the  unsteady  pressure  coefficients 
in  the  vibrstion  of  side  wings  ( NACA64A006  )  . 

1  -  real  part;  2  -  laaginsry  part;  3  -  calculated  by  using 

this  aethod,  4  -  real  part;  5  -  laaginsry  part;  6  -  the 
Tijdeaan  experiaent. 

Pitching  vibration. 

He  Mainly  investigated  the  vibration  situation  of  the 

all-wing  type  for  higher  siaplified  frequencies.  The  results 

l  giver,  in  Fiu  .  3  i u  J 

of  the  real  part  si ag  ree  well.  Ftom  the  results  of  the 
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imaginary  part*  in  Fig.  5(b),  it  can  be  seen  that  the 
relative  errors  of  the  results  in  this  method  compared  to 
the  experimental  results  are  less  than  in  the  LTRAN2  method. 
This  shows  that  when  the  frequency  increases,  the 
iaproveaent  of  the  unsteady  term  is  aainly  reflected  in  the 
iaaginary  parts  which  indicate  the  changes  of  phase  lag  of 
unsteady  pressure,  and  this  will  influence  the  magnitude  of 
the  torque.  The  relation  between  the  tiae-field  pressure 
coefficient  and  the  real  and  iaaginary  parts  of  the 
f requency- f ie Id  is: 

ujI j  +  Io>ca*(u)f)3  (15) 

where  Re  is  the  real  part  and  la  is  the  iaaginary  part. 
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F i i .  5.  The  unsteady  pressure  coefficient  of  the 
NACA64A0 1 0- type  wing  in  pitching  vibration  (x/c=0.25). 

1  -  experiment;  2  -  this  method;  3  -  reference  12};  4  - 
real  part;  5  -  imaginary  part. 

In  this  method,  the  time  step  used  in  the  unsteady- 
flow  calculation  is  usually  2-5  times  larger  than  in 
L  TRAN 2  method  (u£t=6*  -15*).  The  calculation  is  stable  and 
the  region  of  calculation  is  smaller  (same  as  the  above 
network  used  in  steady  flows).  Therefore  the  computer  time 


can  obviously  be  reduced.  The  convergence  standard  of  the 


unsteady  calculation  was  taken  (Atflaax  <  10"*.  Usually  3 
6  circulations  are  required  to  reach  the  standard. 

Our  work  has  obtained  the  sincere  help  of  Professor  Lu 
Shijun  and  Comrade  Li  Jungfu.  The  authors  would  like  to 
thank  thea. 
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Abstract 

A  numerical  method  for  calculation  of  two-dimensional,  inviscid.  steady 
and  unsteady,  small-disturbance  transonic  flows  is  presented.  A  monotone 
■witch  ADI  algorithm  is  used  to  solve  the  governing  equation.  The  a- 
ppropriate  unsteady  terms  are  added  to  the  boundary  conditions  on  body 
surface  and  wake.  An  analytical,  asymptotic  solution  and  the  nonreflecting 
far-field  boundary  conditions  are  chosen  for  steady  and  unsteady  case, 
respectively.  A  nonumiform  mesh  spacing  and  relatively  large  time-step 
are  allowable.  The  developed  computer  program  can  be  applied  for  solving 
both  the  steady  and  unsteady  flows.  Hence,  it  can  be  more  efficient  to 
use. 

The  computational  results  obtained  from  the  present  method  for 
some  examples  have  provided  a  closer  agreement  with  experimental  res¬ 
ults  compared  to  other  methods. 


APPLICATION  OF  DUAL  HOLOGRAM  INTERFEROMETRY  TO  DENSITY  FIELD 


DETERMINATION  IN  WIND  TUNNELS 

Ding  Hanquan,  Chen  Xingyi,  and  Chen  Shouzhi 
Beijing  Institute  of  Aeronautics  and  Astronautics 

1.  INTRODUCTION 

Recently,  pulse  double-exposure  hologram  technology  has 
been  used  in  the  parameter  detection  of  high-speed  wind 
tunnels.  That  is:  two  instant  hologram  pictures  of  the  flow 
field  are  recorded  on  one  negative.  If  there  is  no  relative 
movement,  in  all  the  components  under  the  test,  between  the 
two  exposures,  then  the  interference  pattern  will  reflect 
the  relative  parameters  correctly.  However,  during  the  wind- 
tunnel  test,  because  of  the  strong  disturbance  of 
vibrations,  it  is  very  difficult  to  avoid  disturbance  in 
these  components;  therefore  there  will  be  background  fringes 
in  the  interference  pattern,  and  this  will  lead  to 
deformation.  Furthermore,  in  order  to  measure  and  read  the 
numbers  of  the  interference  fringes  correctly,  the  finite- 
fringe  method  is  usually  used  in  the  recording,  that  is: 
before  the  second  exposure,  one  introduces  a  wedge-like 
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phase  body,  or  the  light  is  deflected  by  a  small  angle,  to 
obtain  the  reference  background  fringes.  However,  it  is 
difficult  to  control  the  orientations  and  distance  between 
the  background  fringes  and  a  sandwiched  holograph  can  mend 
above  shortages  to  a  certain  degree,  but  the  adjustable 
degrees  of  freedom  are  fewer.  We  have  developed  a  negative 
film  shelf  with  six  degrees  of  freedom.  By  using  dual- 
holograph  technology,  the  density  distributions  of  the  axi- 
symmetric  shock  waves  with  the  Mach  number  Ma=2  are 
quantitatively  measured  in  a  real  wind-tunnel  environment, 
with  four  mathematical  models  [1,2].  The  experimental  data 
agrees  well  with  the  AGARD  cone-like  flow-chart  curves. 

2.  THEORY 

Assuming  light  passes  through  the  transverse  cross- 
section  of  the  axi-symmetric  flow  field.  (Fig.  1),  when  the 
flow  field  is  in  a  disturbed  state  (the  refractive  index  is 
n)  and  then  changes  to  a  non-disturbed  state  (the  refractive 
index  is  n„),  the  difference  of  the  refraction  index  is 
related  to  the  displacement,  at  the  point  (y,  z),  of  the 
fringes  of  the  interference  pattern  of  the  double-exposure 
holograph  as  follows: 

S(y,  y,  O -njdx 
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Fig.  1.  A  simplified  diagram  for  calculation 
of  shock-wave  flow  field. 

1-cone;  2-disturbance  zone;  3-shock-wave  lines;  4-shock-wa\e 
surface;  5-cone  surface;  6-light. 
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For  the  cross-section  of  a  given  disturbance  radius  R,  we 


have : 


S(^) 


-n:., 


(1) 


This  is  an  Abel  equation.  The  transformation  of  this 
equation  is: 


■  (  r 


2l  f*  jfOO 
*  J  r  //-r* 


(2) 


Assuming  « (o)  —  b,*-v,  •  From  (1)  and  the  step-function 

method,  we  obtain: 


N-  l 


-tfySCr,)-  2  rJKi'P) 


i  *  i 


B  (  i ,  i  ) 


(  »  -  0,  1,  2,-,AT-l) 


where 


(  3  ) 


According  to  (1 )and  the  linear-function  method,  we  obtain: 

N-  i 

2 

it  -  i  ♦  i 


w 


ysiji)-  2  »o 


V,— 


HUT 


(•-1.2,  3  ,  •••»  N  —  1) 


(4) 


where 


A(  i  )-(  •  +  1  )y/Zi  +  l  -,*ln  ±±  l  +'/2i  ±1 

i 

B(  i ,  >*)-(>*  +  i  +  i  )»-i*-2ii»/]:rri*-+(,i  _  j  )l/(n  _  x y_r^i 

_ _ _ (k+k'h‘Z7*-  )» 

c*  +  »v(hi  )*— i*3c  h  - 1  + ,/fin 
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+»*la 


and 


N-  1 


V»*E ~tfr S(  0  )—  2  2  v» 


i*  «i 


According  to  (2)  and  the  plane-function  method,  we  obtain: 


2  X 


JV-  1 
2 


»» »  • 


2£  +  1 


(5) 


According  to  (1)  and  the  exponential-series-function  method, 
we  obtain: 


1 

2i? 


So  CijL+j  111 

Dm  2‘m 


(6) 


*  «  # 


where  \ — jp-  .  During  the  calculation,  n  is 

taken  up  to  3 . 

According  to  the  Gladstone-Dale  equation,  we  have: 


Pi-P.-VK 


<7) 


For  air  and  ruby  lasers,  the  constant  K=2 . 2582X10_4-m3/kg . 


3.  EXPERIMENT 


All  the  optical  components  and  the  film  shelf  for  the 
holograph  were  attracted  to  a  vibration-isolated  steel-plate 
shelf  by  using  magnetic  bases.  This  makes  all  the  components 
of  the  holographic  recording,  except  the  laser  device, 
connected  as  a  rigid  body. 

Fig.  2  is  the  dual-holograph  shelf.  It  is  a  mechanical 
apparatus,  which  can  hold  two  holograpic  dry  plates,  has 
three  degrees  of  freedom  of  transverse  movement  and  three 
degrees  of  freedom  of  rotation,  and  can  adjust  the  relative 
position  of  the  two  dry  plates  precisely. 


Fig.  2.  The  dual-holographic  shelf. 
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Si-He-Ne  laser  device;  S 2-Ruby-laser  device;  Mi,  M2,  Mj- 
total-ref lection  mirrors;  Li,  Lj-collimating  mirror;  D- 
diffusing  plate;  B-dividing  plate;  m-cone  model;  T-wind- 
tunnel  region;  H-dry  plates  of  holograph. 

The  model  is  a  metal  cone  with  a  half  angle  of  15®.  The 
axis  direction  is  same  as  the  direction  of  the  incoming 
flows;  therefore  an  axi -symmetric  flow  field  is  formed  at 
the  testing  region. 

Fig.  3  is  a  schematic  diagram  of  the  optical-path 
system.  After  the  resonant  cavity  of  S2  and  the  optical  path 
are  adjusted  by  using  Si,  two  dry  plates  of  the  holograph 
are  inserted  to  the  dual-holographic  shelf.  An  exposure  is 
made  before  the  wind  blows  and  the  static  holographic 
pattern  I0  and  II0  are  obtained.  Then  I0  and  II0  are 
detached  and  two  other  dry  plates  are  inserted.  Then  the 
wind  blows  and  another  exposure  is  made.  Thus  the 
holographic  patterns  I  and  II  are  obtained.  By  detaching  I 
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(or  II)  and  combining  it  with  II0  (or  I0) ,  the  dual- 
holographic  pattern  is  obtained  by  pattern  rebuilding.  The 
II  (or  I)  remaining  on  the  dual-holographic  shelf  was 
exposed  to  the  static  field  and  then  the  usual  double- 
exposure  holographic  pattern  is  obtained.  It  is  used  for 
comparison  with  the  dual-holographic  pattern. 


•1 23456 


ftftflUi  ('? 


Fig.  4 

(a)  finite-fringe  interference  pattern;  (b)  a  schematic  diagram 
of  the  measurement  and  reading. 

(l)-the  cross-section  of  reading. 


4.  DATA  MEASUREMENT,  CALCULATION,  AND  RESULTS. 


The  cross-section  for  data  reading  was  taken  at  the 


position  which  is  2/3  cone  length  f roa  the  cone  vertex  (Fig. 
4).  The  length  between  Rc  and  R  were  divided  into  N  segaents 
and  altogether  N+l  ■easuring  points  were  taken.  The  distance 
between  two  neighboring  seasuring  points  is  W=(R  -  Rcl/N,  Rc 
and  R  are  the  radius  of  the  cone  at  the  reading  point  and 
radius  of  the  disturbance,  respectively.  By  using  a 
projector,  we  Measured  Rc  =  0.7100  cm,  R  =  1.7750  cat .  With  N 
=  20,  The  serial  nuabers  i  at  Rc  and  R  are  15  and  35, 
respect i ve 1 y . 


r  <  ft. 


Fig.  5.  A  coMparison  of  the  experiMental  results  and  the 


AGARD  curve. 


Isualiy,  we  always  define  the  displacement  of  the 
interference  fringes  at  the  intersection  point  of  the 
reading  cross-section  and  Bhoek  waves  as  S=0  for  the 
calculation.  As  a  matter  of  fact,  this  is  true  only  at  the 
first  exposure  to  the  static  field  which  is  an  undisturbed 
cons t ant - f 1 ow  field,  (That  is;  wind  blows  without  the  model, 
the  dens i t >  of  the  incoming  flow  is  ,  and  the  refractive 
index  is  n„).  When  the  static  field  is  still  (the  densitv 
and  refractive  index  of  the  environmental  air  are  i> 0  and  n0, 
res pe<  tivelyi;  the  S  value  at  that  po i n t  should  be  ( n , - 
n  i  l.  /  -  ,  (1  is  the  transverse  distance  of  the  testing 

rej i  n  .  \s  we  have  defined  S=0  for  that  po int  ,  therefore, 

*•  i  t  h*  - 1  ’ho  s’atii  field  is  still  or  moving;  we  can  treat  it 
as  moving,  thus  n0  and  e0  in  the  above  formulas  should  be 
changed  into  n,  and  P  *  . 

The  interferes  e  fringes  are  coded  according  to  Fig. 

4 ' h i .  Wh  »-n  the  background  fringes  at  the  undisturbed  region 
are  paral  lei  to  the  ■  •  r os s - se<  ■  t  ion  of  the  reading,  the  S  1  , 

>f«  at  the  intersecting  points  of  each  fringe  of  the 

distur bed  region  and  this  cross-section  will  be  the  code 
values.  If  there  is  a  slope  angle  between  the  background 
fringes  and  the  cross - »ec t i on ,  the  measured  S’  should  be 
corrected  in  order  to  obtain  the  real  value  S. 

As  a  results  of  the  measurement  ,  the  total  pressure 

2  7 


r 


p0:2  atm,  the  total  temperature  T0=281  K.  By 
equation,  r0;.  ~  5  1  kg/m3  .  When  Ma«=2,  by  using 
•>,/P0  =  0-2301  .  Then  P«=0.5772  kg/m3. 

The  density  near  the  inside  flank  sides 


us i ng  the  state 
the  table, 

l  s  : 


Pa 


(  h  +1  )/(  k  -  1  ) 
-  1  +  2  /(  k  -  1 


Substitute  the  ratio  of  heat  capacities  k  =  1  .  4  ,  Ma«=2, 
and  shoik-wave  angle  ^=33. 7*  into  above  formula,  the 
ini’rement  at  the  shot  k-wave  surface  is: 

Pa-P./P.- AP/P.-0.19 

The  experimental  data  are  shown  in  Fig.  5.  It  can  be 
seer,  that  the  experimental  results  agree  well  with  the  AGARI' 
curve.  Especially  the  s t ep- f unc t l on  method  (with  simple 
calculations  and  small  accumulated  errors)  agrees  quite 
well.  This  shows  that  the  theory  and  method  proposed  in  this 
paper  is  appropriate  for  investigating  the  density 
distribution  of  ax i - s ymme t r l c  high-speed  air  flows. 
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This  paper  'presents  the  make  use  of  the  dual  hologram  interferometry 
for  obtaining  finite  fringes  holographic  interferogram  of  the  a  x  isvmmetric 


flow  field  with  shock  waves,  and  the 
four  mathematical  models. 


density  distributions  computed  by 
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